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Abstract 

The method of partial Fourier transform is used to find explicit for- 
mulas for two remarkable fundamental solutions for a generalized Tri- 
comi operator. These fundamental solutions reflect clearly the mixed 
type of the Tricomi operator. In proving these results, we establish 
explicit formulas for Fourier transforms of some functions involving 
Bessel functions. 



1 Introduction 

Consider the generalized Tricomi operator 
(1.1) P = yA+ ' " 



dy 2 ' 



"Partially supported by CNPq (Brazil) 
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where A = Y] n , 7—^ is the Laplace operator. Our aim is to find, via partial 

3 ~ ox 1 * 

Fourier transform with respect to ), fundamental solutions 

relative to an arbitrary point (a, 0) on the hyperplane y = in R n+1 , that is, 
distributions that are solutions to 

Pu = 8(x — a,y), 

where 5(x — a,y) is the Dirac measure concentrated at (a,0). Since P is 
invariant under translations parallel to that hyperplane, it suffices to consider 
the case when the Dirac measure is concentrated at the origin. 
If n = 1, then ( |1.1| ) is the classical Tricomi operator, 

Q2 Q2 

(1-2) ^ = te + 



dx 2 dy 2 

It is known that, for this operator, the equation 9x 2 + 4y 3 = defines the 
two characteristic curves that originate at the origin. They divide the plane 
in two disjoint regions, namely, 

D+ = {(x, y) eM 2 : 9x 2 + % 3 > 0}, 

the region "outside" the characteristics and 

D- = {(x, y) e M 2 : 9x 2 + Ay 3 < 0}, 

the region "inside" the characteristics. Note that _D_ is entirely contained in 
the hyperbolic region y < 0. 

In the paper 0] it was shown the existence of the following two funda- 
mental solutions relative to the origin: 

:i . 3) F + {x,y) = \ 3 ' 22/3 * 1/2r PV3) 



and 

(1.4) F_(x,y) 



elsewhere 
3r(4/3) 



\9x 2 + 4y 3 |- 1/6 inD_ 



2 2 /3 7r i/2r( 5 / 6 ) 
elsewhere. 



2 



The support of F_, the closure of D_, is, except for the origin, entirely 
contained in the hyperbolic region (y < 0), while the support of F + , the 
closure of D + , consists of the whole elliptic region (y > 0), the parabolic 
region (x-axis), and extends to the hyperbolic region up to and included the 
characteristic curves. The method used in || to prove these results was based 
upon the property that solutions to the equation Tu = are invariant under 
the dilation d t (x,y) = (t 3 x,t 2 y) in R 2 . 

In [0] we proved formulas (|1.3| ) and (|1.4|) by using partial Fourier trans- 
form with respect to the x variable. In this paper we extend the results of Jl| 
to the generalized Tricomi operator ( |1.1| ). Since the dimension of the space 
variable x is now n > 1, technical difficulties in evaluating Fourier transforms 
involving Bessel functions do occur. We circumvent them by calculating in- 
tegrals of the type 



I e (a,b) 



r x JJat)J v {bt)dt 



(Section |). As a result, we obtain for the operator ( |1 . 1| ) the following fun- 
damental solutions: 



jn-2 



■p/n I s 

12 sJ (9|x| 2 + 4y 3 )i"t i n D r , 



F+{x,y) 



22/3^/2 r( 2 /3) 



elsewhere, 



whose support is the closure of the region D\ = {(x, y) € R™ +1 : 9|x| 2 +4?/ 3 > 
0} and 



F-(x,y) 



3 r( Q ) i i i n o.ln „ 

7TtI 9 M +% 3 |3"2 i n D n _ 



2 2 / 3 W 2 r(| - f ) 




elsewhere 



supported by the closure of the region Z?" = {(x, y) £ M n+1 : 9|x| 2 -|-4?/ 3 < 0}. 

These fundamental solutions clearly generalize formulas (|1.3|) and ( [L.4|) . 
In addition, we also show that no matter the parity of the dimension n, 
F_(x,y) does not have support on the boundary of the region D_ in the 
hyperbolic region and hence, in contrast with what happens for strictly hy- 
perbolic operators - like the wave operator - the Huyghens principle does 
not hold for the generalized Tricomi operator. 
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It is well known || that, for the wave operator, the regularity of the 
fundamental solutions degenerates as the dimension n increases: a locally 
constant function, when n = 1; an absolutely continuous measure relative 
to the Lebesgue measure, when n = 2; a measure carried by the surface of 
the forward light-cone, when n = 3; and so on. However, for the generalized 
Tricomi operator they always remain locally integrable. 

The plan of this paper is as follows. In Section || we reduce, via partial 
Fourier transform, the original problem to an equivalent one of finding fun- 
damental solutions for a second order ordinary differential equation and show 
how these can be represented in terms of Airy functions or Bessel functions. 
In Section |3| we obtain explicit formulas for Fourier transforms of the func- 
tions Id^ld), lei^(lel), and lel^lel), £ = (6, • • • , £n), n > I. These 
formulas are used to construct the fundamental solution supported by D'l 
(Section f|) and the fundamental solution supported by (Section In 
the Appendix the reader will find the definitions of the Bessel functions J u (z), 
I„(z), K v {z), and N u (z), and a list of properties that are needed throughout 
this work. We also recall the definition of hypergeometric functions and some 
of their main properties. 

The method of finding solutions, via partial Fourier transform, to the 
equation 

(1.5) Pu = yAu+ — = f 

with, say, / G C^°(SL n+1 ) is a natural one although not particularly new. In 
the monograph H, R. J. P. Groothuizen exhibits a solution u in terms of 
a Fourier integral operator whose symbol is obtained by the partial Fourier 
transform analysis employed in this paper. He also obtains fundamental 
solutions given by Fourier integral operators. Our results are more precise 
since we give explicit formulas for the fundamental solutions and these are 
tempered distributions. Consequently, we can obtain more transparent rep- 
resentations for the solution u to the equation ( |1.5[ ) as a convolution of / 
with any of the fundamental solutions here described. 

2 Preliminaries 

Consider the more general problem of finding fundamental solutions for the 
operator ( |1 . 1| ) relative to an arbitrary point (0, b) on the y-axis. That is, one 
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wishes to find distributional solutions to the equation 

(2.1) yAF + — F = 6{x) <g> 6{y - b). 

Partial Fourier transform with respect to x reduces this problem into finding 
fundamental solutions to the ordinary differential equation 

(2.2) Fyy-yltfF = 6(y-b), 

where F denotes the partial Fourier transform of F, |£| 2 = ^™ =1 £j, and 
S(y — b) is the Dirac measure concentrated at b. 

One way of solving (|2.2j ) (see fH) consists in selecting two linearly inde- 
pendent solutions Ui(£,y) and U2(£,y) to the homogeneous equation 

(2.3) u yy -y\i\ 2 u = Q 

so that their Wronskian at y = b (which in the case under consideration is 
the same as the Wronskian at y — 0) is equal to —1 and in defining 



(2.4) Ht,y;b) 



U 2 {i,b)U x {^y) \iy>b 
Ui(£,b)U 2 {Z,y) if y<b. 



It is a matter of verification that F(£,y;b) is a solution to ( |2.2|) . 

Equivalently [ffl, one can obtain a fundamental solution to the equation 
( |2.2j ) by finding two linearly independent solutions to the homogeneous equa- 
tion ( p.3|) : Ui(£,y; b), defined for y > b, and v,2(£,y; b), defined for y < b, so 
that 

(i) the limit of u\ as y —>■ b+ equals the limit of u 2 as y — > b— 

and 

(ii) the limit of u^ y as y — > b+ minus the limit of u 2t y as y — > b— is equal 
to 1. The function F(£, y; b) is now defined by 



(2.5) F(£,y;b) 



ui(£,y,b) if y > b 



In what follows we will use at our convenience either one of these two expre- 
sions for F(£, y \ b). 
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Returning to the original problem (|2.1|), if we can choose U%(£,y) and 
U2(£,y), in formula Q2.4Q , or ui(£,y;b) and w 2 (£, &), in formula ( |2.5|) , so 
that F(£,y;b) is a tempered distribution with respect to £ = (£i,---,£n)> 
then its inverse Fourier transform of -F(£, ?/; 6) defines a fundamental solution 
to the operator 

We now proceed to find linearly independent solutions to the homoge- 
neous equation (|2.3| ). The change of variables z = |£| 2//3 ?/ transforms (|2.3|) 
into the classical Airy's equation u" — zu = 0. Two linearly independent 
solutions to that equation are Ai(z) and Bi(z) respectively called the Airy 
functions of the first and second kinds. It is known that these two func- 
tions can be represented in terms of Bessel functions of order ±1/3 as follows. 
If | argz| < (2?r/3), then 



. Ai(z) = ^[/_ 1/3 (f, 3 / 2 ) - I 1/3 (F /2 )] = i(i ) 1/2 K 1/3 (lz^), 

(2.6) 

' 5z(,) = (|)V2[ / _ 1/3 (| 2 3/ 2)+/i/3( | 2 3/ 2)] . 

If | &rgz\ < (2-7T/3), then 



(2.7) 



Ai(-z) 



1/2 



[J_ V 3(f2 3/2 ) + J V 3(f<Z 3/2 )], 



5z(-^) = (|)V2[J_ 1/3 (|,3/2 ) _ Ji/3( | z 3/ 2)] . 

In the Appendix I, where a brief review of Bessel functions is presented, 
we show that the following relations hold: 

3-2/3 3-4/3 
(2-8) Az(0) = — — — , Ai'(0) -- 



and 
(2.9) 



Bi(0) 



r(2/3) : 

3-I/6 



Bi'(Q) 



r(4/3) 

3-5/6 



r(2/3)' v y r(4/3)' 

As a consequence, the Wronskian of Ai(z) and Bi(z) evaluated at z = is 
(2.10) W(^,fli(s))|^ = l/7r. 

Indeed, we have 

3-2/3 3-I/6 



W(Ai(z),Bi(z)) u=0 



r(2/3) r(2/3) 

3-4/3 3-5/6 



r(4/3) r(4/3) 



2 ■ 3- 3 / 2 _ 1 

r(2/3)r(4/3) ~ tt : 
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because 

(2.11) r(2/3)r(4/3) " 



33/2 • 

We now choose the following two linearly independent solutions to the 
equation 



-7r|e|" 2/3 ^(|e| 2/3 &)^(le| 2/3 y) if V > b 
-n\i\-^Ai{\i\^b)Bi{\^y) ily<b 



U 1 (^y) = ^\- 1/3 At(\e /3 y) and U 2 & y) = -yft\Z\- 1/3 Bi(\Z\ 2 ' 3 y), 

and note that, by virtue of fl2.10| ), the Wronskian of U\(£,y) and U2(£,y) is 
equal to —1. Next, according to (|2.4j), the distribution 



(2.12) F(Z,y;b) 



is a fundamental solution to the ordinary differential equation ( |2.2| ). A fun- 
damental solution F(x,y;b) to the generalized Tricomi operator (11.1]) and 
relative to the point (0, b) would then be the inverse Fourier transform of 
F(^,y; b) whenever that Fourier transform exists. 

We do not know how to obtain an explicit formula (or formulas) for the 
inverse Fourier transform of F(£, y; b) when b 7^ 0, a problem that merits to 
be investigated. We conjecture that when b < 0, that is, the point (0, b) is in 
the hyperbolic region, there exists two fundamental solutions that converge, 
as b — > 0, to the fundamental solutions F + (x,y) and F_(x,y) defined, re- 
spectively, by the formulas ( |5.6| ) and ( [4.2|) . The two fundamental solutions 
described in the monograph || do not seem to satisfy these requirements. 

However, when b = 0, we will show in the following sections how to obtain 



from formula (|2 . 1 2j ) , as well as formulas similar to it, explicit expressions for 
fundamental solution to ( |1.1|) . 

The case n = 1. For sake of completeness and in order to motivate our 
work in the forthcoming sections we briefly present the results of the paper 



Hfl, for the Tricomi operator (|1.2j ), that is the case when n — 1. If b — and 
we take into account the values of Ai(0) and Bi(0) as given by ( J2.8|) and 
Q2.9Q , then formula (|2.12|) simplifies as follows: 



(2.13) F((,y) 



-It I -2/3 

s4(273) Bi «| 2/3 ^ lf ^°' 
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where, for simplicity, we wrote F(£, y) for F(£, y; 0). Now the inverse Fourier 
transform of both expressions on the right-hand side of ( [2.13|) can be explic- 
itly calculated. To see this, first consider the representations of Ai(z) and 
Bi(z) in terms of Bessel functions as given by formulas ( |2.6|) and ( |2.7| ): 

1/2 



and 



Bim 2/3 y) 



Ai{\i\ 2 "y) 



\i\ 2/ \-y) 



1 

7T 

1/2 



K| 2/3 y 



K 



[«/-i/3(|iei(-2/) 3/2 ) - ^/3(|iei(-?/) 3/2 )]- 



Next by introducing the change of variables s = (2/3)y 3 ^ 2 , whenever y > 0, 
and t = (2/3)(— y) 3 ^ 2 , whenever y < 0, rewrite formula Q2.13|) as 

:^) 1/3 k 1/3 (s\z\) if y >o 



a: 



(2.14) F(£,3/)=< 



where a and j3 are constants given by 

1 



(2.15) 



a 



and 



7T 



2 1 /3 3 i/3r( 2 /3) — ^ 2 1 /3 3 5/6 r ( 2 /3)' 
Theorem 3.1 of [|IJ proves that F(x,y), the inverse Fourier transform of 
F(£,y), is then 



If. 

2 



a linear combination of the two fundamental solutions F + and F_ respec- 



tively defined by formulas Q1.3Q and (|1.4|) . To prove this theorem, one re- 
lies on known formulas (see 0,0]) for the inverse Fourier transforms of the 
functions l^ 1 / 3 J 1/3 (|£|), |£|- 1/3 J-i/ 3 (|£|), and |£|- 1/3 ^i/ 3 (|£l), formulas that 
need be generalized to the case n > 1. 



3 Fourier transforms involving Bessel func- 
tions 

As we have indicated in the previous section, we are going to need ex- 
plicit formulas for the inverse Fourier transforms of the following functions: 
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lei^oei), ler^cei), and ler^ciei), where e = (6, . . . , e„), « > i- since, 

when n = 1, explicit formulas for the Fourier transforms are known and can 
be found, for example, in || and J7|, we concentrate on the case n > 1. 

1. The inverse Fourier transform of l^-fGd^l). 
Assume that |Re v\ < 1/2 and define 

(3.i) j-MiewKDioo = 7^ / e^)|er^(iei)de- 

The assumption on the real part of i/ secures convergence of the integral at the 
origin. On the other hand, in view of the asymptotic behavior of i^Q^D for 
large values of |£|, the integral converges absolutely. By introducing spherical 
coordinates, rewrite the integral on the right-hand of ( |3.1|) as 

oo p 

r n +u - lK ^ r ^ e i{rx ' u) duj} dr. 



Since n 

(3.2) / du = p^-Jn x {r\x\), 
Js n -! \rx\ 2 

it follows that 

| T |l-f roc 

(3.3) ^- l m u KM\W) = [ 77-^ r r ^Jn_MA)K»{r)dr 

(27T) 2 J 



We now quote the following result found in Watson's treatise |1 I 
Lemma 3.1. If Re(fi + 1) > \Re v\ and Re b > \Im a\ then 



POO 

(3.4) / V L+u+1 J^at)K u (bt)dt 

Jo 



(2aY(2b) v T(^ + v + l) 



(a 2 + b 2 Y+ u+1 
From this lemma, we obtain the following result 
Theorem 3.1. If \Re v\ < 1/2 then 

(3.5) r-vmtm*) = — (! + NT 

7T2 
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Proof. The integral in ( |3.3| ) is the same as the integral in ( |3.4| ), where 
/i = n/2 — 1, a = \x\, and 6=1. Clearly the conditions of the lemma are 
satisfied and Theorem |3.1| follows at once. □ 
We remark that, when n = 1, (13.51) becomes 



7T2 

a formula found in ||] and . 

2. Inverse Fourier transforms of an d ICl^^dCl)- 

As before, assume that |Re z/| < 1/2. Formally, the inverse Fourier trans- 
form of \£\"J V (\Z\) is 

^ 1 [ier^(iei)](x) = -^- / e^iewieode. 

In general, the integral diverges at oo and so we introduce a converging 
factor e~ e ^' and take a limit as e — > 0. Since is locally integrable, 

it defines a tempered distribution and the limit exists in <S'(IR n ), the space 
of tempered distributions on R n . Thus the precise meaning of the inverse 
Fourier transform is 

(3.6) r-^wuizMx) = / e^-^ier^dei)^. 

[ATT) e-HJ J Rn 

By introducing spherical coordinates, one can see that the integral on the 
right-hand can be written as 

-er n+v-l j ( \ r / J{rx,u) 



e- er r n+u - L J,(r){ / e l{rx ^ duo} dr. 



Since 



(3.7) / e ^^=-^ T Jn„ 1 (r|x|), 



\rx\ 2 



it follows that 



I II— ^ roc 

(3.8) F-'m'MmW = lim / e^r? +"Jn _ x (r|x|) J„(r) dr. 

(2tt)2 ^oJ 2 
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In a similar manner we also have 



(3.9) ^[ItrJMim = r e-^r^J^ MA)Ur)dr. 

Both integrals appearing in formulas ( [3.8|) and ( |3.9|) are particular cases 
of the following integral 

POO 

(3.10) I e (a,b)= e~ et r x J^at)J„(bt) dt, 

Jo 

where a and b are positive real numbers, A, /i, and u, complex numbers such 
that Re(/x + v + 1) > Re(A). This integral is a variant of the discontinuous 
integral of Weber and Schafheitlin studied by Watson in his treatise on Bessel 
functions 



Lemma 3.2. As e — > I e (a,b) tends, in the sense of distributions, to either 
one of the following limits: 



(I 



if < b < a, or 

,0^ a^r(^±l) » + y _A + l ^-A-i/ + l a* 

v- 1 ') 2 ^+ 1 r(// + i)r(to±i) 1 2 ' 2 '^ +i ' &2 ^ 

z/ < a < 6. 

Proof. We adapt to our situation the proof of the Weber-Schafheitlin theo- 
rem as found in section 13.4 of Watson's treatise It consists of expanding 
the integrand in ( |3.10 ) in power series of b, and passing to the limit as e — > 0. 

1. Consider the case when < b < a. If we replace b by z, the integral 
( |3.10| ) is an analytic function of z when Rez > and |Imz| < e. Introduce 
new constants a, /3, and 7 defined by 

2a = /i + i/-A + l, 2/? = z/-A-// + l, 7 = z/ + l 

or, equivalently, 

A = 7 — a — /3, [i = a — j3, u = j — 1, 
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and rewrite ( |3.10 ) as follows: 

POO 

(3.13) I e (a,z)= / e" e 't a+/3 ~ 7 J Q _ /3 (at)J 7 _ 1 (^)rft 



-et T , ^fsr- (-i) m (z/2y +2m -H a+ P +2m - 1 1 

e J Q _^(at)^ > — — ■ V 

l ' m!l it + to) J 



m=0 

m— 1 ^oo 



m!T(7 + to) 



The interchange betwen the integration and summation signs is justified be- 
cause, when \z\ < e, the series 

f^ n ™!r( 7 + to) Jo 



is absolutely convergent. 

2. We now evaluate the last integral on the right-hand side of (|3.13|) . By 
expanding J a _p(at) in power series and integrating term by term we obtain 



POO 

/ e- rf J a _^(at)t a+/3+2m - 1 dt 
Jo 

00 ( 1 \ k (n io\ a -P+ 2k r°° 

k=o k\r(a-p + k + i)J 

(-l) fc (a/2)^+ 2fc r(2a + 2m + 2A;) 

" h\V(rv — /-? —I— jEr* —I— T ^2a+2m+2fc ' 



fc=0 



By using the duplication formula T{2z) = 2 2z ~ l ix~ 1 l 2 Y(z)T( y z + |) together 
with the notation {z) k = z(z + 1) • • • {z + k — 1) = r(z + A;) /T(z) we rewrite 
the last expression as: 



e~ et J^^at)^^ 2 " 1 - 1 dt 



2\ k 



(a/2) a ^ T(2a + 2m) ^ (a + m) k (a + to + \) k 
(e 2 ) a + m T(a-P + 1)^ k\(a-P + l) k V c 2 

k=0 

(a/2) a ~P T(2« + 2to) 1 a 2 , 

' F(a + m } a + m + -;a-/3+ 1; — -). 



(e 2 ) Q+m r(a-/3+ 1) v 1 ' 1 1 2' ^ ' c 2 
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By using formula ( 6.1 6| ) in the Appendix II: 



F(a, b; c; z) = (1 - z) a F(a, c-b;c; -), 



we obtain that 
(3.14) 



-ct 



o+/3+2m-l 



dt 







F(a + m,a + m H — ; a — (3 + 1; 



a 2 + e< 



(a 2 + e 2 ) a+m T(a- (3 + 1 
Substituting ( [04] ) into ( |3TT3D we get 
(3.15) / e (a,z) = 



^ -l m (z/2 7+2m ^ 1 (a/2 a ~ / T 2a + 2m _ 1 fl 1 

> v , , wo — ^rf^ — ^ ^-F(a+m,a+m+-;a-3+l; 

^ m\T 7 + m a 2 + e 2 "+ m r (ar-/3+l) V 2' P ' 



m=0 

whenever |z| < e. 

3. Following [11], one can show that (|3.15 ) is valid provided that z satisfies 
the conditions 



a 2 + 



Re(z) > 0, |Im(z)|<e, 



< \/a 2 + e 2 - e. 



Let 5 > be small enough so that 6 < a 2 + S 2 — S, and take < e < S so 
that we also have 6 < Va 2 + e 2 — £■ In ( p,15| ) we may now let z = b and, when 
this is done, one can show, by the method of majorants, that the resulting 
series converges uniformly with respect to e, < e < 5, and therefore, as 
e — > 0, the limit of the series is equal to its value at e = 0. Thus 



(3.16) 

oo 

= £ 

m=0 



lim I e (a, b) 



1 m (6/2 T+ 2 ™- 1 (a/2 a -^r 2a + 2m ) „, 1 

; ,)/ ; r \! ' K - r — -F a+m, a+m+-; a- 0+1: 1 . 

m!r 7 + m a 2Q+2 "T a - (3 + 1 v 2' M ' ; 



By formula ( p. 15 ) in the Appendix II, we have 
F(a + m, (3 — m; a — (3 + 1; 1) 



r(a-/3+l)r(l/2) 



r(l-/3-m)r(a + m + i)' 
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On the other hand, 



r(l/2)r(2a + 2m) = 2 



2a+2m-l 



r(a + m)T(a + m + -) 



T(l-P-m) 



•l) m 7TCSc(7T/3) 



T((3 + m) 
T(1-(3)T((3) = 7TCSC(7T/?), 



therefore 



r(2a + 2m)r(l/2) 



2 2a+2m - 1 T(a + m)T{f3 + m) 
r( 7 + m)r(l -(3- m)T(a + m + |) ~ (-1)"T(1 - /3)r(/3)r( 7 + m) 

Substituting these formulas into Q3.16D , we get 



(3.17) 



lim I e (a, b) 



E 

m=0 



feT- 1 r(a + m)r(,5 + m) 



2 7-a-/3 a a+/3 m ir(i - /3)r(/3)r(7 + m) 



2 \ m 



2T- a -^a a +^r(7)r(i - /3) 



F(a,(3;j; 



Finally, returning to the constants /z, z/, and A, we obtain the expression 
( 3.11 ) in the first part of the lemma. 

4. In the case < a < b, we proceed in an analogous manner to obtain 
the expression (|3.12|) in the second part of the lemma. □ 

We now use Lemma B.2I to evaluate the inverse Fourier transforms of 



If r^Gf |) and \£.\~ V J»(\£\) respectively defined by formulas (p|) and (p|). 

Theorem 3.2. The inverse Fourier transform of \£\ V J V (\£\) is the distribu- 
tion defined by 



(3.i8) r-^trumw 



sin 



mr 2*T(§ + u) 



7T2 



+ 1 



.r 



— sin (un 



< Id < 1, 



1 < \x\ 



7T2 
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Proof. The integral on the right-hand side of ( |3.8| ) corresponds to the inte- 
gral in ( ft.lOQ where 

n n 
/i = — — 1, v — v, A = — — — z^, a — \x\, b—1. 

If 1 < \x\, then from formula ( |3.11| ) we obtain for the limit in ( |3.8|) : 

e~ er r 2 +u Jn_ 1 (^ r \x\) J v {r) dr = 

F(- + z/,z/ + l;z/ + l; 



2~^\x\^ +l Y(u + l)Y(-u) v 2 
From the known relation for hypergeometric series 
(3.19) F(a, b; c; z) = (1 - z) c ~ a ' b F(c -a,c-b;c;z 

we have that 

ii | 
F(- + i/, i/ + 1; i/ + 1; 



x 



\x 


2 


- 1 




X 


2 



2 ' ' ' |af' 
On the other hand 

i> + i)r(-i/) = 

Thus, for 1 < \x\, the above limit is equal to 

e - er ri +v Jn^{r\x\) J v {r) dr = 

= -Sin(7TZ/) — ^ (|x| 

Z 2 7T 

If < | a; | < 1, then formula (|3.12|) yields 

POO 

(3.21) lim / e- T ri + "j!i_ 1 (r\x\)J„(r)dr 



1)" 



2 

) 



z|3-T(2 + i/ 



n n n 



— Ff u — ■ — ■ It 

2-f-«T(|)r(l - I) l 2 ' 2' 2' 1 



- 2 



sin 



?M±!0 (1 _ wr .-.. 



2 7 2-2-^vr 
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after making the replacements 



r£)r(i-^ 

K 2 J v 2' 



7T 



and 



. n 



n n 



F(- + u, . . 
v 2 2 2 



sin(wr/2) 



'1 — Id 



The last expression results from ( |3.19| ). Substituting ( |3.2(J| ) and (3.21) into 
(|3.8|) we obtain (|3.18 ) and the theorem is proved □ 

By reasoning in the same manner the following result, whose proof is left 
to the reader, holds. 

Theorem 3.3. The inverse Fourier transform of \£,\~ u Ju{\£,\) is the distri- 
bution defined by 



(1 - \x 



(3.22) r-^trumi* 



2 v ix^T{y-\ + 1) 



0. 



< Id < 1, 



1 < \x\. 



Remark. From theorem [3.3| we immediately derive a formula for the 
inverse Fourier transform of |£| v </-i/(|£|), namely, 



1 — Id 



(3.23) ^m-j-Mmix 



2-"7rir(l - v - 1; 



< \x\ < 1, 



0. 



1 < Id 



Since 



7T 



r(i-,--)r(, + -) = - ^ 



we may rewrite Q3.23D as 

(3.24) r\\tfJ- v (\Z\)\{x) = 

" 2^sin[(z/+f)vr]r(z/+f) 



7T2 



+ 1 



X 



\2\-u- 



2 , o < Id < 1, 



0, 1< \x\. 

3. The inverse Fourier transform of ^"-/V^d^l). 
Theorems 3^2 and 3^3 imply the following result 
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Theorem 3.4. The inverse Fourier transform of \£\ V N V (\£\) is the distribu- 
tion defined by 



(3.25) ^-MKWieDK*) 



,n7r,2T(f +v) 



cos(— ) n +1 (1- MT^ 0<|d,<l, 



7T 2 



— COS Z/7T 



2T(§ + z/) 

2.4-1 



(>| 2 - lj-y-", 1 < \x\. 



Proof. Recall that 



J,(|e|)cos(^7r)- J_„(|£|) 



Thus, from formulas ( |3.18| ) and ( |3.24| ) we get 

cos(z/7r) 2 ! T(^ + z/) 

■ \ 5Ti (1-H) 2 ~ 



Sm(^7T] 7T2 



+ 1 



-cob(i/7t) \ 2 + | J (N 2 -l^ 



< Id < 1, 



1 < Id 



7T2 



Sin (U7T) 



-+1 



Now 



0. 



7T Ti7T 

sin[(z/ + 77)] 71 " — sin(z/7r) cos(— ) 



1 < Id 



,nir. 



COS Z/7T sin 



and so formula ( |3.25| ) follows at once. □ 



4 A fundamental solution with support in D^l 



We return to the problem (|2.1|) and, in order to obtain a formula for a fun- 
damental solution with support in the hyperbolic region, we have to modify 
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formula ( |2.13| ). Consider the function 

r 3 2 / 3 r(4/3) 
(4.1) F.( { ,„)= 2VS 



for y > 0, 



where t = (2/3)(— y) 3 ^ 2 . From formula (|6.10|) in the Appendix I, it follows 
that the limit of y) as y — > 0— is equal to zero and so condition (i) in 
Section ^| holds. Also it follows from the same formula that 

^{(^) 1/3 ^l/3(^l)}| y =0 ' '"' 



^y -w\suj\y=*> 3 2 /3r(4/3)' 

hence, the y-derivative of y)aty — 0— is equal to —1, and so condition 
(ii) is also satisfied. By calculating the inverse Fourier transform of F_(£,y) 
we have the following result: 

Theorem 4.1. The inverse Fourier transform of F_(£,y) is the distribution 

|9|a;| 2 + % 3 |3"2 in 



(4.2) F_(x,y) — < *T(i"S) 

[ elsewhere, 

where = G M n+1 : 9|a;| 2 + % 3 < 0}. It zs a fundamental solution 

for the operator whose support is the closure of the region .D™. 

Proof. We must evaluate the inverse Fourier transform of (t / \^\Y^ 3 Ji/ 3 (t\^\) . 
Recall that if G(x) is the inverse Fourier transform of /(£), then (l/a n )G(x/a) 
is the inverse Fourier transform of /(a£). By appying formula (|3.22| ) for v = 
1/3, we obtain 



^ 1 [(Vlel) 1/ V 1/3 (t|el)]( 



x 



i 



ta n ( \x 



2 1 /3 7r f r(| - 2) 



i n 



3 2 

, < \x\ < i, 



0, t < \x\ 
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2 1 /3 vr im 

\3 2 



— (t 2 - |x| 2 ) 3 2 , < \x\ < t, 



0, £ < \x\. 

3 n_ § i_« 

— |9|x| 2 + % 3 | 3 2 in Z?" 



2V3 7r f r(| - f ) 

0, elsewhere. 

Multiplication by the constant 3 2/3 r(4/3)/2 1/3 in formula flOD yields 
(|4.2j ) which proves the theorem. □ 

Remarks 1. If n — 1, then F_(x,y) coincides with the distribution 
defined by formula ([[21). 

2. For all values of n the support of F_(x,y) is the closure of the region 
D™. This follows from results in [|J about the generalized function where 
P(rc) is the quadratic polynomial 

p( \ — 2 , _|_ 2 _ 2 _ _ 2 

p, q > 1. See also the example in about the distribution Pf. s x , where s 
denotes the hyperbolic distance in ~R N , N > 1. 



5 A fundamental solution with support in D\ 

Consider the function F + (£,y) defined by 

[ 7-(^) 1/3 ^i/ 3 (^|) if2/>0 

(5.1) F + (Z,y) = { f 

{ ^-(^) 1/3 iV-i/3(^l) if 2/ < 0, 

where s = (2/3)y 3//2 , t = (2/3)(— y) 3 ^ 2 , and the constants 7 and S are respec- 
tively given by 

2 2 / 3 2tt 

(5.2) 7 = — — — - - and 5 



3 4 / 3 r(2/3) 2 1 /3 3 4/3 r ( 2 /3)' 

and where AL1/3 is the Neumann function defined by ( |6.4j ). By using the 
formulas in the Appendix I, it is a matter of verification that the conditions 
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(i) and (ii) in Section ^ are satisfied. Thus its inverse Fourier transform, 
denoted by F$(x,y), defines a fundamental solution for the operator ( [Lip . 
The following theorem gives an explicit expression for F*{x,y). 

Theorem 5.1. The inverse Fourier transform of F + (£,y) is the distribution 
defined by 



(5.3) F\x,y) 



3 w ~ 2 r(f-|; 

22/3^/2 r(2/3) 



in D\ 



.wr. 2^3 
— cos(— ) 



n-2 pfn _ I % 



2 7 W 2 T(2/3) 



|9|x| 2 + 4y 3 | 3 2 i n ^ 



Proof. 1. We start by evaluating the inverse Fourier transform of (s/|^|) 1//3 ii'i/3(s|^| 
s = 2w 3 / 2 /3, y > 0. Since K v {z) = K_ u (z), it follows from formula ( p.5Q that 



^(«/iei) 1/s JM«fci)] = ^| |} (* 2 + M 2 ) 



1 _ n 
3 2 



7T 



After multiplying by the constant 7 in ( |5.2j ) and reintroducing the variables 
x and y, we can see that the right-hand side of the last expression is equal to 



(5.4) 



yn— 2-p/n 1' 



2 2 / 3 W 2 r(2/3) 



_L _ n 

(9|x| 2 + % 3 ) 3 2 , y >0, 



that is F tt (x, y) in £>!» n {y > 0}. 

2. Next we evaluate the inverse Fourier transform of (t /\£\y/ 3 N^i/ 3 (t\£\) , 
t = 2(— y) 3 / 2 /3, y < 0. From formula ( [3.25] ) we obtain 



cos 



,n7r 2-V3r(g-i; 



(t 2 - |x| 2 )3 2, < |x| < t, 



i\ 2- x /3 r (| _ i 



7T 2 



^(|0;| 2 » 



■ > 1 » 

t Z 5-7. 



£ < Ixl. 



After multiplying both sides by the constant 8 in ( |5.2P and reverting to the 
variables x and y, we can see that the right-hand side of the last expression 
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can be written as 



(5.5) 



n— 2r(n 1' 



gn-iwn 



1 n 



2 2 / 3 W 2 r(2/3) 



2 — ^(9|a;| 2 + V) 3 2 , inD™n{y<0} 



^Tr^^n-ar (n_I) , a 
cos — / 2 ; . 3J 9d 2 + 4y 3 3 2 , inF 

v 2 ; W 2 r 2/3 111 y 1 



and so (|5.5|) coincides with F$(x,y) in Z?" D {?/ < 0} U ZP. Therefore, from 
(|5.4j) and (|5.5| ) we obtain ( |5.3|) and the theorem is proved. □ 

Remarks 1. In view of the results in [Q relative to the generalized 
function (see Remark 2 after the proof of Theorem |4.1| ) the supports of 
the distributions 



(9|x| 2 + % 3 )3"2 an d \9\x\ 2 + 4y 3 



i 



that appear in formula ( |5.3| ) are the closures of Z)™ and Z)™, respectively. In 
other words, for no value of n can these supports be just the boundaries of 
these regions. 

2. If n is odd, then cos(rar/2) = 0, and we rewrite formula (|5.3| ), using 
the notation F + (x,y) instead of F^(x,y), as follows: 



>n-2 V( n _ I"\ 

5 ^2 3^ /nU|2 i ,<„,3< 1 " 



(5.6) F + {x,y) 



XP+V)3-2 in/}™ 



2 2 / 3 W 2 r(2/3) 

elsewhere. 



This is a fundamental solution whose support is D™. In particular, if n = 1 
we obtain formula (|1.3|) . 

3. If n is even, then F"(x, y) is not necessarily identically zero in Z)™ and 
its support may be the whole of M n+1 . Suppose that n = 2k, k > 0, and let 
us compare the constant in formula ( |5.3|) , relative to the region D™, to the 
constant in the expression ( f4.2|) of F~(x,y). Let 

(_ 1 )fc+i 2 i/3 3 2(fc-i) r(k - ±) 
a = - — — 



7T 



r(2/3) 



be the constant in (|5.3|) and let 
B 



3 2fc T(4/3) {-\) k+1 Z 7k - 1 Y{k-\) 



2 2 / 3 7r fc r(| — fc) 2 2 /% fc r(2/3) 
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be the constant in ( |4.2j ). Since 3 A — 2B = it follows that the distribution 



(5.7) F+(x, y) = 3F«(x, y) - 2F_(x, y) 

is now a fundamental solution for the operator ( |1.1[ ) supported in and we 
have for this F + (x, y) the same expression as that of ( |5.6| ). In conclusion, for 
all values on n we always get two fundamental solutions: one whose support 
is the closure of D 1 } and another whose support is the closure of D™. 



6 Appendix 



I. Bessel functions 

The function J u (z) of a complex variable z defined by 

(6.1) Uz) = Y, 2 v+2r rT ( l/ + r + 1 y l*l<°°> |arg^|<7r, 

is called the Bessel function of the first kind of order v. 

We also need the Bessel functions I„{z) and K„(z) defined by 

00 z u+2r 

(6.2) F(z) = — — — ; — ; -, \z\ < oo, I argd < tt, 

and 

7T CSC I 7/7T ) 

(6.3) K u {z) = ±-J-{I_ v ( z ) - I v ( z )l ^0,±1,±2,... 

as well as the Neumann function 

J v {z) COs(z/7r) — J- U (z) 



(6.4) N v {z) 



sm w 



Note that throughout this work, we only deal with Bessel functions of 
order ±1/3. Recall that Ai(z) was defined in formula ( |2.6|) by 



Ai(z) = ^[/-!/3(^ 3/2 ) -/l/3(^ 3/2 )] = ^(|) 1/2 ^/3(^ 3/2 ) 

If we set s = |z 3//2 , then we may rewrite Ai(z) as 

(6-5) Ai(z) = ^ S 1/3 [Li/ 3 ( 5 ) - h/3(s)]. 
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From the series expansion of I v {z) it follows that 

(6-6) s^i. 1/3 (s) = 2 _ 1/3 r (2/3) + ^57^) + • • • 

and 

s 2/3 ^8/3 

(6 ' 7) S 3/l/3(s) = 2V3r(4/3) + 2V3r(7/3) + " ' 

Consequently from ( p.5[ ), fl6.6|), and ( |6.7|) we obtain 

3-2/3 



Ai(0) 



r(2/3)' 



the first expression in formula (2Ji). Similarly, by differentiating Ai(z) and 
setting z = 0, we get the second expression in ( |2.8| ) 

3-4/3 

Ai'(0) - 



r(4/3)- 

In an analogous way, recall that Bi(z) in Q2.6| ) was defined by 

^) = (|) 1/2 [/-l/3(^ 3/2 )+/l/3(^ 3/2 )]. 

If we set s = f z 3 / 2 , then the last expression becomes 



3 

I 

2V331/6 ' 



(6-8) Bi(z) = T^m^l 1 -^ 8 ) + h/s(s)}, 



and again from (|6.6|) and (|6.7|) we obtain 



3-1/6 3-5/6 
Bi{0) = ^r--, Bi'{0) " 



r(2/3)' v ' r(4/3)' 



which are the two expressions in ( |2.9| ). 

Finally, from these results we obtain the value of the Wronskian of Ai(z) 
and Bi(z) at z = 0: 

2 • 3~ 3 / 2 

W(Ai(z), Bi(z))\ . = ,,,,,, = 1/vr, 
v v h v ;;u=o r(2/3)r(4/3) 7 

because r(2/3)r(4/3) = 2vr • 3- 3 / 2 . 
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For future reference, we also need the following expansions derived from 
the series that defines J v (z): 

1 t 2 

( \ 1/3 T ( \ 

(6,9) f 1/3(t) = 2-V3r(2/3) ~ 25/3r(5/3) + • • • 

and 

f 2/3 f 8/3 

II. Hypergeometric series and functions 

Let a, b, and c be arbitrary complex numbers and let z be a complex 
variable. The power series 



71 

) 



„. , , (a, n)(b, n) z 

6.11 Fa, 6; c; * = V 1 7 A ^ 

c, n) n! 

n=0 v ' 7 

where 

(6.12) (a, 0) = 1, (a, n) = v . , = a(a + 1) • ■ • (a + n - 1), n = 1, 2, . . . 

T(a) 

and we assume that c 7^ 0, —1, —2, . . . , is called a hypergeometric series. It 
is known |12[ that the series (|6.11|) is a solution, valid near z = 0, of the 
hypergeometric equation 

d?u du 

(6.13) z(l - z)— + {c - (a + 6 + l)z\— -abu = 0, 

dz 2 dz 

for which every point is an ordinary point, except 0, 1, and 00, that are 
regular singular points. 

If either a or 6 is a negative integer, then the series ( |6.11 ) terminates; 



if c is a negative integer, the series is meaningless because all terms after 
the (1 — c) th have a zero denominator. As it is known |J, it is possible to 
redefine the series (|6 . 1 1|) so that it still is a solution of the hypergeometric 
equation. We exclude this possibility from our considerations because, in the 
cases that interest us, c is never a negative integer. 

The hypergeometric series is absolutely convergent for \z\ < 1 and so 
defines, in the open disk, an analytic function of z which is regular at z — 0. 
The point z — 1 is however a branch point and if a cut is made from 1 
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to +00 along the it can be shown ||12j| that series can be continued 

analytically and defines an analytic function throughout the cut plane that 
we still denote by F(x, b; c; z). If Re(c) > Re(6) > 0, this analytic extension 
can be represented by Euler's formula 

(6.14) F(a, b; c; z) = T{b ^_ fc) J ^(1 - b) M (l - tz)~ a dt, 

for I arg(l — z)\ < 71. 

In general, the hypergeometric series ( |6.11|) diverges for \z\ = 1. However, 
if Re(c — a — b) > 0, we have absolute convergence for \z\ = 1. Moreover, 

(6.15) F( a ,6;c;l) = ^ r ( C - a -g . 

r(c — a)T(c — 0) 

Hypergeometric functions satisfy among themselves quite a number of 
important relations of which we just list the following two that we used in 
Section [|: 

(6.16) F(a, b- c; z) = (1 - z)- a F(a, c-b;c; —?—) 

z — 1 

and 

(6.17) F(a, b; c; z) = (1 - z) c_ °- 6 F(c - a, c - 6; c; 2). 

For a complete list of such relations, the reader should consult Erdely H]. 
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